Instability and control of a periodically-driven Bose-Einstein condensate 
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We investigate the dynamics of a Bose-Einstein condensate held in an optical lattice under the 
influence of a strong periodic driving potential. Studying the mean-field version of the Bose-Hubbard 
model reveals that the condensate becomes highly unstable when the effective intersite tunneling 
becomes negative. We further show how controlling the sign of the tunneling can be used as a 
powerful tool to manage the dispersion of an atomic wavepacket, and thus to create a pulsed atomic 
soliton laser. 
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Introduction - The spectacular experimental progress 
in confining Bose-Einstein condensates in optical lattice 
potentials has provided a powerful tool for investigating 
many-body quantum dynamics. Such optical potentials 
are extremely clean and controllable, and together with 
their long decoherence times, this allows the observation 
of many coherent lattice phenomena which are extremely 
challenging to study in other solid-state systems. One 
such effect is "coherent destruction of tunneling" (CDT) 
observed very recently in atomic systems 0, H], in 
which a periodic driving field acts to renormalise the 
tunneling between lattice sites. This control over the 
dynamics of the condensate is achieved without altering 
any of the parameters of the optical lattice, and has nat- 
ural applications to quantum information, since it pre- 
serves the system's coherence. However, it is crucial to 
know the stability of the condensate during its time evo- 
lution, particularly the presence of dynamical instability, 
in which deviations from a steady state grow exponen- 
tially with time. The case of a static potential was ana- 
lyzed in Ref.jl], and studied experimentally in Refs.[l,@], 
and it was found that dynamical instability occurs above 
a certain critical quasimomentum. Later work Q ex- 
amined the case of a uniformly accelerated lattice, and 
found the non-intuitive result that dynamical instability 
was enhanced in the limit of low acceleration. 

In this paper we build upon this approach to analyze 
the richer and more complex case where the optical po- 
tential is periodically rocked. To achieve this we carry 
out the stability analysis about the Floquet states ^J] of 
the system, which are the appropriate generalization of 
energy eigenstates to the case of a time-periodic Hamil- 
tonian. An important point for experiment is to mini- 
mize or avoid instabilities, and so we first find the crit- 
ical interaction strength at which dynamical instability 
occurs. We then connect this with the behavior of the 
effective tunneling J^g, and show how manipulating J^g 
via CDT can be used to control the dynamics of the con- 
densate, providing control over matter-wave dispersion 
and thus allowing the creation of bright solitons. 

Method - A system of cold bosons held in an optical 
lattice can be described very accurately by the Bose- 
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Hubbard Hamiltonian 
■f^BH = - J ^ (aL^n + H.c. 

{m,n) m 

(1) 

where al^/am are the boson creation/annihilation oper- 
ators, and Um = al^Om is the number operator. The 
properties of the system are governed by the hopping 
parameter J, and the Hubbard interaction U which de- 
scribes the potential energy between two bosons occupy- 
ing the same lattice site. An extremely valuable means 
of studying and controlling such systems is to accelerate 
the lattice by varying the phase-difference between the 
two laser beams forming the standing wave potential. In 
the rest frame of the lattice this acceleration manifests 
itself as an inertial force which effectively "tilts" the po- 
tential. If instead of a uniform acceleration the lattice is 
periodically accelerated and decelerated, it is possible to 
produce a potential that oscillates periodically in time. 
Hi = K cosiot'^^ rnUm- where K and co parametrise its 
amplitude and frequency respectively. 

In order to study the stability of the driven conden- 
sate we will first pass to a mean-field description of this 
model, analogous to the Gross-Pitaevskii equation, and 
then linearize about the ground-state to obtain the Bo- 
goliubov equations for the condensate excitations. We 
first write the Heisenberg equations of motion for the 
boson operators a„, and then take the classical field ap- 
proximation and treat them simply as c- numbers a„. It 
is then straightforward to show that the classical ampli- 
tudes obey the equation of motion 

i— ^ = -J(a„+1 +Q;„_i)+g|Q:„| + X cos cjt nQ„ (2) 
ot 

where for convenience we have scaled the interaction as 
g — U/N. Note that we also take Ti — 1, and will measure 
all energies in units of J. 

To simplify the analysis we use periodic boundary con- 
ditions. In the limit of large lattice sizes, however, the 
choice of boundary conditions does not affect the under- 
lying physics, and we will later use Dirichlet boundary 
conditions to simulate the time-evolution of the conden- 
sate. In the absence of interactions, the eigenstates of the 
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system will simply be plane waves a„ = exp [inp] where p 
is the angular momentum. With this in mind, we take as 
a trial solution = exp [i {n(j) + 9)], where (p and are 
functions to be determined. Substituting this solution in 
EqEl yields the result 

K 

4>{t) = p sintji (3) 

e{t) = 2 J (cosp S{t) ~ sinp C{t)) ~ gt (4) 

where the functions S{t)/C{t) are defined in terms of 
Bessel functions as 
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(6) 



Since the Hamiltonian of the system is periodic in time, 
the Floquet theorem dictates that the solutions of the 
time-dependent Schrodinger equation can be written in 
the form exp [lei] u(£), where e is termed the quasienergy, 
and u[t) is a T-periodic function called the Floquet 
state. To obtain the quasienergies we thus simply have 
to extract the terms from the solution which are not T- 
periodic, giving the result 



e(p) = 2 cosp J Jo {Kjuj) + g. 



(7) 



In the absence of the driving the quasienergies thus form 
a normal single-particle bandstructure, the interaction g 
acting merely to shift the entire spectrum. The driv- 
ing then acts to renormalise the width of the spectrum 
by the Bessel function j Tp, as was previously observed 
in a theoretical analysis [lOj of semiconductor superlat- 
tice systems. In Fig|TJi we show numerical results for the 
quasienergies of an 8-site system, obtained directly from 
the time-evolution of the system, which beautifully cor- 
roborate the expected behavior. In particular, when the 
Bessel function becomes zero, the spectrum collapses to 
a point and the system will manifest CDT. 

In order to analyze the dynamical stability of the 
ground state (p = 0), we now introduce a perturbation 
Oinit) = a" (t)(l + exp [i(7n] exp \-~iqn}\), where 
is the unperturbed solution, and q and w are the 
momentum and energy of the excitation. We then lin- 
earize Eql2] about this solution to obtain the Bogoliubov 
de Gennes equations for u(t) and v(i) 



dt 



u{t) 
v[t) 



C{q,t) 



u(t) 
v{t) 



(8) 



where the elements of the matrix C{q^ t) are given by 

= 4Jsin ((7/2) sin (g/2 — iC/tj sinwt) -I- g, 
^i2{q,t) = 5 = ~^2i{q,t), 

^22{q,t) = — 4Jsin (q/2) sin (g/2 + if/tjsintjt) — 5. 



We can note that, similarly to the Hamiltonian, the oper- 
ator C{q,t) is T-periodic. Consequently we can also ap- 
ply the Floquet theorem to describe the time-evolution 
of the quasiparticle excitation (w, u). To find the corre- 
sponding Floquet states, we numerically evolve Eq[5]over 
one period of driving, using the 2x2 identity matrix as 
the initial state. The result of this procedure is the single- 
period propagator U. The eigenstates of U are then the 
excitation Floquet states, while its eigenvalues are related 
to the excitation quasienergies via = exp [— jTe^]. 

The symmetries of combined with the normalization 
condition obeyed by the quasiparticle excitation, |wp — 
I up = 1, allow the characteristic equation to be written 
in the particularly simple form Q 



A± = Rc [Un] ± [Un] - 1. 



(9) 



The condition of dynamical stability can now be recast 
in terms of the eigenvalues of C{q, t), in analogy with the 
use of Lyapunov exponents in classical mechanics. The 
solution is stable against an excitation with momentum 
q if there is no quasienergy with a positive imaginary 
component, or equivalently, if Re [Un] < 1. This suggests 
a simple scheme to map out the stability zones of the 
driven system. For a given choice of K and oj we select a 
value of g and scan over the range of q. If Re [J/n] < 1 for 
all values of q we can declare that the system is stable for 
these parameters, and that to induce instability we need 
to increase g to a higher value. In this way a standard 
bisection scheme can be used to locate the instability 
boundary, gc- 

Results - We show the results of this procedure 
for two driving frequencies in FiglTjD. We can first 
observe that for K/uj — it is possible to di- 
rectly diagonalise C{q) to obtain the result e± = 

±2V2 sin {q/2) ^ 2 J2 sin^ {q/2) + Jg. As expected, this 
duplicates the familiar result for the Bogoliubov excita- 
tions of an undriven, stationary condensate ll| . It is also 
clear from this expression that dynamical instability will 
not occur if the interaction is repulsive, since the eigen- 
frequencies will not become complex unless the product 
(J g) < 0. Accordingly, as K/uj ^ we can see from 
Figdb that the value of gc diverges. For larger K/uj the 
value of gc then rapidly drops, passing through a broad 
local minimum before again diverging as K/uj approaches 
2.4048 - the first zero of Jq. This corresponds to the on- 
set of CDT; as the effective tunneling is reduced, the 
dynamics of the condensate is suppressed, and stability 
is regained. 

Passing through the zero of Jo, we can see that dy- 
namical stability is then abruptly lost. In this region the 
condensate becomes dynamically unstable for any pos- 
itive value of the interaction. We can obtain some in- 
sight into this effect from Eql?! by defining an effective 
tunnehng, Jgfj = JJq{K/u)). When K/uj is increased 
from 2.4048, the Bessel function changes sign and Jgff 
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FIG. 1: (a) Quasienergy spectrum of the driven mean-field 
Bose-Hubbard model, for an 8-site system with uj — 16. The 
width of the spectrum is modulated by the Bessel function 
J'o{K/uj) and displaced by the interaction energy g — 0.5, 
in full agreement with the analytical solution Eq[7] (b) Plot 
of the critical interaction, g^, at which the system becomes 
dynamically unstable, for driving frequencies uj — 4 and 8. 
This quantity diverges at K/u = and at the Bessel function 
zeros (vertical dashed lines), and is zero when Jgg is negative, 
for which dynamical instability occurs for any positive value 
of g. 

becomes negative. The physical significance of this sign- 
change has been observed previously in experiment 
where it caused the momentum distribution function to 
be discretely shifted by tt. This occurs because the tun- 
neling, as well as being renormalized in amplitude, ac- 
quires a phase- factor of exp[i7r] [T^- Accordingly, if we 
view the driving field as acting simply to renormalise the 
tunneling, the product ( Jgff g) now becomes negative and 
dynamical instability can indeed occur. A similar effect 
would occur if instead g were made negative (for example, 
by using a Feshbach resonance), for which the condensate 
would become attractive and thus unstable toward col- 
lapse. We can note that in Ref. [1] the condensate was 
close to non-interacting, and thus avoided this instabil- 
ity. When Jefj becomes zero at K/uj = 5.52 we can see 
that gc again diverges due to CDT. The same pattern of 
behavior then repeats. Fig[T)3 also shows that gc scales 
quite accurately as gc ^ w^. Surprisingly, the zone of 
stability thus becomes wider at high driving frequencies, 
although the acceleration of the lattice is much larger. A 
similar feature was seen in the analysis of an accelerated 
condensate [t!], which found an increasing propensity to 
dynamical instability in the limit of low acceleration. 



We have so far used an ideal flat optical lattice po- 
tential. In experiment, however, an additional harmonic 
trap potential is usually present which can substantially 
modify the dynamics of the system, and introduce new 
effects. To investigate this, we now apply an additional 



quadratic potential V = krf, where ri is measured from 
the center of the system. We initialize the system in the 
ground state of the the mean-field Hamiltonian ^ in 
the presence of the trap, and then displace the trapping 
potential by a distance of 25 lattice spacings, thereby 
exciting the condensate into motion. 

In Figl2^, we show the evolution of the condensate 
in the absence of the periodic driving. The condensate 
makes a periodic oscillation of constant amplitude, very 
similar to the center of mass motion observed in exper- 
iment Q. The period of oscillation is governed by the 
intersite tunneling, or equivalently, by the condensate's 
effective mass (m* cx J~g), and thus allows these quan- 
tities to be measured directly. In Figllb the system is 
subjected to a driving with K/uj — 2. Clearly the oscilla- 
tion period has increased, corresponding to the expected 
reduction of the tunneling by j7o, which can alternatively 
be interpreted as an enhancement of the effective mass. 
Increasing K/oj further to the first zero of j7o produces 
CDT, and so the system remains frozen in its initial state 
(Figlll:). In this case the effective mass has become infi- 
nite. 

A further increase of K/u) means that the Bessel func- 
tion changes sign, and thus Jeff becomes negative. We 
have seen already that this sign-change has a significant 
effect on the dynamical stability of the ground-state, and 
as we show in Figl2li it has an equally dramatic effect on 
the dynamics of the condensate. Instead of oscillating, 
the condensate now rapidly accelerates away from the 
center of the trap. The reason for this becomes evident 
when we examine the terms of the Hamiltonian. Making 
Jeff negative is clearly equivalent to time-reversed evolu- 
tion with a positive tunneling, but with reversed signs for 
the trapping potential and non-linearity g. Thus when 
Jeff changes sign, the condensate behaves as an attractive 
condensate in an inverted potential, and so is quickly ex- 
pelled from the center of the trap. 



Controlling the parameters of the periodic driving field 
thus allows us to tune the effective mass of the conden- 
sate to be positive, negative, or infinite. A particularly 
important application of this control is the production of 
solitons. Solitons remain stable by balancing dispersion 
with the interparticle interaction, and thus bright soliton 
solutions of the Gross-Pitaevskii equation demand either 
using an attractive interaction, or the use of complicated 
staggered phase-imprinting techniques. Making the effec- 
tive mass negative, however, would allow bright solitons 
to be created in repulsive condensates, without having 
to reverse the sign of the interaction or requiring phase 
imprinting. 

To investigate this possibility, we revisit a proposal 
made by Carr and Brand in Ref. T3!| for manipulating a 
trapped condensate to produce a train of solitons. The 
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FIG. 2: Dynamics of a condensate in a parabolic trap. At 
t — the trap is shifted by 25 lattice spacings to induce the 
condensate into motion. System parameters are g — 0.1 and 
OJ = 4. (a) In the absence of driving, the condensate sloshes 
periodicaUy from side to side of the trap with a well-defined 
frequency, (b) For K/u — 2.2, J^ff is reduced in amplitude, 
and so the frequency of the oscillation is correspondingly re- 
duced, (c) For K/u = 2.40, the tunneling is completely sup- 
pressed since Jgg ~ 0, and the time-evolution of the conden- 
sate is thus frozen, (d) For K/ui = 3, Jgg becomes negative. 
The trap potential thus appears inverted to the condensate, 
which is rapidly expelled from the center. 

procedure consists of two parts; the scattering length of 
the condensate is changed from positive to negative while 
at the same time the trap is inverted to become expul- 
sive. As we have seen previously, both of these processes 
can be accomplished simultaneously by tuning the sign of 
Jgff from positive to negative. To verify this, the trapped 
system was again initialized in its ground-state, and the 
amplitude of the driving slowly increased from zero to 
K/u! = 2.0 so that the system adiabatically followed to 
the ground-state of the renormalized Hamiltonian. As 
can be seen in the top panel of Fig|31 the resulting density 
distribution is gaussian in form. The amplitude of the 
driving was then abruptly altered to K/uj = 3.8, revers- 
ing the sign of J^g . The resulting time evolution of the 
condensate shown in Fig[3] clearly displays the transfor- 
mation of the condensate into a train of solitonic pulses, 
with the solitons forming first at the edges of the con- 
densate. 



Conclusions - In summary, we have shown how an os- 
cillating driving potential can be used to renormalise the 
effective tunneling, Jgg, of a Bose-Einstein condensate. 
This parameter crucially determines when the conden- 
sate is dynamically unstable, and when Jgfj is negative 
the condensate becomes unstable for any positive value 
of interaction. This is particularly important for the de- 
sign of experiments, since in such parameter ranges the 
interaction must be set to an appropriately small value. 



We have also shown how manipulating J^fl- in this way 
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FIG. 3: Time-evolution of a condensate held in a parabolic 
trap, subjected by a periodic driving field. System parameters 
are g = O.I, and lj = 4. At t = the amplitude of the tunnel- 
ing is suddenly changed from K/uj = 2.0 to 3.80, to reverse 
the sign of Jgg. (a) At t = the condensate has a gaussian 
density-profile, (b) t — 200T, the condensate has spread, and 
has developed oscillations at its edges, (c) t — 400T the oscil- 
lations have deepened into localized pulses, similar to those 
seen in Ref . [l3l] . and the initial state has almost completely 
converted into a soliton train. 



may also be used as a novel tool to to control the dy- 
namics of a condensate, in a complementary way to the 
well-known method of controlling the interaction via Fes- 
hbach resonances. This both extends the possibility of 
manipulating the condensate to systems which do not 
possess convenient resonances, and provides a new means 
to investigate the interplay between nonlinearity and dis- 
persion, notably the production of solitons. 

The author was supported by a Ramon y Cajal Fel- 
lowship. 
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